iT邦幫忙

2018 iT 邦幫忙鐵人賽
DAY 25
0
AI & Machine Learning

玩轉資料與機器學習-以自然語言處理為例系列 第 25

分群演算法實作-區辨不同演算法的意義

  • 分享至 

  • xImage
  •  

前言

這篇文章將承襲上一篇文章,視覺化各個演算法的成果,並附註說明,讓大家更能夠了解,每個演算法運作的原理,及其用處。以下將用五個資料集,進行實測: flower, noisy_circles, noisy_moons, blobs, gaussian_quantiles。這原本是用來做分類的,所以紅色跟藍色就是他的類別,然後這個資料集為了方便視覺化,所以只設計出兩個維度。

Import

import matplotlib.pyplot as plt
from planar_utils import plot_decision_boundary, sigmoid, load_planar_dataset, load_extra_datasets
import numpy as np
import os

from sklearn.cluster import KMeans
from sklearn.cluster import AgglomerativeClustering
from sklearn.cluster import DBSCAN
from collections import Counter

資料準備

X, Y = load_planar_dataset()
name = 'flower'
X = X.T
Y = Y[0]
plt.scatter(X[:, 0], X[:, 1], c=Y , s=40, cmap=plt.cm.Spectral);
plt.title(name+'_original')
plt.savefig(os.path.join('pic', name+'_original'))
plt.show()

noisy_circles, noisy_moons, blobs, gaussian_quantiles, no_structure = load_extra_datasets()
datasets = {"noisy_circles": noisy_circles,
            "noisy_moons": noisy_moons,
            "blobs": blobs,
            "gaussian_quantiles": gaussian_quantiles}

datas = [(name, X, Y), ]
for name, dataset in datasets.items():
    X, Y = datasets[idx]
    datas.append((name, X, Y))
    
    print(X.shape)
    print(Y.shape)

    # Visualize the data
    plt.scatter(X[:, 0], X[:, 1], c=Y, s=40, cmap=plt.cm.Spectral)
    plt.title(name+'_original')
    plt.savefig(os.path.join('pic', name+'_original'))
    plt.show()

https://ithelp.ithome.com.tw/upload/images/20171227/201075760AVP7net45.jpg

Flower

這張圖中,你們可以發現kmenas與hierarchical演算法的效果其實差不多,其實都是找比較相近的點,kmeans用比較有效率的方式找,hierarchical則是一群一群慢慢整併。dbscan則透過密度去尋找群組,可以設定eps的效果就不是特別好,因為dbscan比較常用來尋找極端值,而flower這張圖,很明顯不是一整群比較密集配上其他零星的小點,因此,密度比較高的自己就會被分成一群。
https://ithelp.ithome.com.tw/upload/images/20171227/20107576xSyxubuDqY.jpg

noisy_circles

這張圖來看,你就會發現dbscan的效果特別好,因為圓形內部的密度很高,而外部的點密度比較低,自然就會把密度高的分成一群。
https://ithelp.ithome.com.tw/upload/images/20171227/20107576j81E0ivZUA.jpg

noisy_moons

https://ithelp.ithome.com.tw/upload/images/20171227/201075762bbWVSQ7eS.jpg

blobs

https://ithelp.ithome.com.tw/upload/images/20171227/20107576tDnzpsmQaG.jpg

gaussian_quantiles

https://ithelp.ithome.com.tw/upload/images/20171227/201075764I4028VTfG.jpg

不同eps的scan

大家會發現,在eps==0.2的時候,他的分群效果可以很明顯。如果要找極端值,eps設在0.3或0.4會比較剛好。
https://ithelp.ithome.com.tw/upload/images/20171227/20107576uU2sAp3DcU.jpg

https://ithelp.ithome.com.tw/upload/images/20171227/20107576Y4wDOtK6hf.jpg

其他比較圖片的生產

# kmeans
for name, X, Y in datas:
    kmeans = KMeans(n_clusters=n_clusters, random_state=0).fit_predict(X)
    plt.scatter(X[:, 0], X[:, 1], c=kmeans, s=40, cmap=plt.cm.Spectral);
    plt.title(name+'_kmeans')
    plt.savefig(os.path.join('pic', name+'_kmeans'))
    plt.show()
    
# hierarchical clustering: single link
linkage = ['ward', 'average', 'complete']
for name, X, Y in datas:
    hierarchy = AgglomerativeClustering(linkage=linkage[0], n_clusters=n_clusters).fit_predict(X)
    plt.scatter(X[:, 0], X[:, 1], c=hierarchy, s=40, cmap=plt.cm.Spectral);
    plt.title(name+'_single')
    plt.savefig(os.path.join('pic', name+'_single'))
    plt.show()
    
# hierarchical clustering: average link
for name, X, Y in datas:
    n_clusters = len(set(Y))
    hierarchy = AgglomerativeClustering(linkage=linkage[1], n_clusters=n_clusters).fit_predict(X)
    plt.scatter(X[:, 0], X[:, 1], c=hierarchy, s=40, cmap=plt.cm.Spectral)
    plt.title(name+'_average')
    plt.savefig(os.path.join('pic', name+'_average'))
    plt.show()
    
# hierarchical clustering: complete link
for name, X, Y in datas:
    n_clusters = len(set(Y))
    hierarchy = AgglomerativeClustering(linkage=linkage[2], n_clusters=n_clusters).fit_predict(X)
    plt.scatter(X[:, 0], X[:, 1], c=hierarchy, s=40, cmap=plt.cm.Spectral);
    plt.title(name+'_complete')
    plt.savefig(os.path.join('pic', name+'_complete'))
    plt.show()
    
# dbscan
for name, X, Y in datas:
    n_clusters = len(set(Y))
    dbscan = DBSCAN(eps=0.3).fit_predict(X)
    plt.scatter(X[:, 0], X[:, 1], c=dbscan, s=40, cmap=plt.cm.Spectral)
    plt.title(name+'_dbscan')
    plt.savefig(os.path.join('pic', name+'_dbscan'))
    plt.show()

code在這裡


上一篇
分群演算法理論
下一篇
產品標籤分群實作-Word2Vec
系列文
玩轉資料與機器學習-以自然語言處理為例31
圖片
  直播研討會
圖片
{{ item.channelVendor }} {{ item.webinarstarted }} |
{{ formatDate(item.duration) }}
直播中

1 則留言

0
s790502ss
iT邦新手 3 級 ‧ 2022-01-12 19:14:36

請問版主 planar_utils 是要額外安裝的獨立套件嗎?

我要留言

立即登入留言